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ABSTRACT 

Th e hypervelocity sta r SDSS J090745. 0+024507 in the halo of the Milky Way 
galaxy ijBrown et alj f2005) most likely o riginated fr om the breakup of a binary star 
system by the central black hole, SgrA* l)Hillslll988j) . We examine the fate of former 
binary companions to similar hypervelocity stars (HVSs) by simulating 600 differ- 
ent binary orbits around SgrA* with a direct N-body integration code. For some 
orbital parameters, the binary breakup process leads to HVSs with ejection veloci- 
ties that are almost an order of magnitude larger than the velocity observed for SDSS 
J090745. 0+024507. The former companion stars populate highly eccentric orbits which 
resemble the observed orbits for some of the stars nearest to SgrA*. 

Key words: black hole physics-Galaxy:center-Galaxy:kinematics and dynamics- 
stellar dynamics 



1 INTRODUCTION 

Recently, the first hypervelocity star (HVS), SDSS 
J090745. 0+024507, was discovered in the Galactic halo 
{Brown et al J 1200.4 iFuentes et al.ll2005h . This HVS is lo- 
cated at a heliocentric distance of ~ 71 kpc and has ra- 
dial velocity 853 ± 12 kms -1 . Its velocity is over twice that 
needed to e scape the gravitational pull of the Milky Way. 
iHillj il988f) was the first to recognize that a HVS might 
result from a close encounter between a tightly bound bi- 
nary st ar system and the blac k hole at the Galactic center, 
SgrA*. lYu fc Tremaind <|2003T) refined Hills' argument and 
added that HVSs might also be produced by three-body in- 
teractions between a star and a binary black hole. Because 
the existence of a second blac k hole in the Galactic center 
(Hansen & Milosavlievic |2003h is only a hypothetical possi- 
bilitv (ISchodel et al.ll200aT we focus our discussion on the 
disruption of a tightly bound binary by a single supermas- 
sive black hole (SMBH). 

The Keplerian orbits of massive stars within 10 2 -10 4 AU 
from the Galactic center provide strong evidence for the ex- 
istence of a central SMBH with mass ~ 4 x 10 6 M W (e.g. 
[ Ghez et al] l2004 iReid fc Brunthalerl 120041: ISchodel et all 
2003). Since binaries are common in other star forming 
environments, it is only natural to explore the interaction 
between SgrA* and nearby binaries. Both lYu fc Tremaind 
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(120031) and lGould fc" Quillenl 12003) describe instances where 
the tidal disruption of a binary by the SMBH leads to one 
star being ejected into close orbit around the black hole. 
Thus, it is only natural to ask: Is it possible that some of 
these stars are former companions of HVSs? 

In §2 we describe the N-body code and simulation pa- 
rameters that were adopted. In §3 we discuss our numerical 
results for the origin of the HVSs, and in §4 we compare 
the calculated orbits of the bound companion stars to the 
observed stellar orbits near the Galactic center. Our goal is 
not to cover the entire phase space of possible binary orbits 
but rather to examine whether some of the highly eccen- 
tric orbits of observed stars near SgrA* could have resulted 
with a reasonable probability from the disruption of a stellar 
binary. 



2 COMPUTATIONAL METHOD 

In our study we have used the N-body code N BODY writ- 
ten bylAarsetbl lll999l) . and presented in lBinnev fc Tremaind 
(1987). We have tested NBODY0 against later versions of 
Aarseth's N-body codes (such as "triple"), and found the 
results to be identical to within the required precision. We 
adopted a small value of 10 -8 for the accuracy parameter rj, 
which determines the integration step through the relation 
dt = y/ nF / (d 2 F / dt 2 ) where dt is the timestep and F is the 
force. The softening parameter, eps2, which is used to create 
the softened point-mass potential, was set to zero. We treat 
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the stars as point particles and ignore tidal and general rela- 
tivistic effects on their orbits, since these effects are small at 
the distance (~ 10AU) where the binary is tidally disrupted 
by the SMBH. 

We have set the mass of the SMBH to M = 4 x 1O 6 M 
and the mass of each star to m = 3Mq, comparable to 
the estimated mass of SDSS J090745.0+024507 (Fuentes 
et al. 2005). All runs start with the center of the circu- 
lar binary located 2000 AU (= 10~ 2 pc) away from the 
SMBH along the positive y-axis. This distance is compa- 
rable to the inn er scale of the observed distribution of stars 
around SgrA* jEckart fc Genzelll997l: ISchodel et aT1l2003t 
iGhez et al] 12003) 7 allowing the remaining star to populate 
this region after the ejection of its companion. This radius is 
also much larger than the binary size or the distance of clos- 
est approach necessary to obtain the relevant ejection veloc- 
ity of HVSs, making the simulated orbits nearly parabolic. 
We used the same initial distance for all runs to make the 
comparison among them easier to interpret as we varied the 
distance of closest approach to the SMBH or the relative 
positions of the two stars within the binary. 

We chose initial binary separations of a = 0.05 or 0.1 AU 
because they provide ejection velocities in the range of inter- 
est 1 for the above parameters. Sig nificantly wider binarie s 
would give lower ejection velocities dGualandris et al 
Much tighter binaries would not be easily disrupted by the 
black hole, or may coalesce to make a single star before inter- 
acting with the SMBH. The size of a main sequence star of 
a few solar masses is ~ 0.01AU, and so binaries tighter than 
~ 0.02AU are precluded because the two stars will develop 
a common envelope and eventually coalesce. 

In the Galactic disk, about half of all stars form in bina- 
ries o r small multiple systems (see e.g. iDuauennov fc Mavorl 
Il99ll) . with roughly equal probability per loga rithmic in- 
terval of separations. dP/ dln(a) = const (e.g. lAbtHl983l 
lHeacoxllT99a: iLarsonlbool . In the Galactic center environ- 
ment, the maximum binary separation is limited by the tidal 
force of SgrA* at the distance d where the binary is formed 
(fo r conditions that enable sta r formation near the SMBH, 
sec Milosa vlievic fc Loebll2004T) . Since the mass of the black 
hole is ~ 10 6 times larger than that of a star, this implies a 
maximum binary separation less than (10 _6 ) 1//3 = 10~ 2 of 
the initial distance d. For d = 2 x 10 3 AU, the upper limit 
on the binary separation would be 20AU (or smaller if the 
tidal restriction applies during the formation process of the 
binary). If we assume a constant probability per ln(a) for 
0.02 < a < 20AU, then the probability of finding a binary 
in the range of a = 0.05-0. 1AU is substantial, ~ 10%. 

We have found that the initial phase of the binary orbit 
plays a crucial role in the outcome. Therefore, we sampled 
cases with initial phase values of: -75, -60, -45, -30, -15, 0, 15, 
30, 45, 60, 75, and 90 degrees. As initial conditions, we gave 
the binary system no radial velocity but a tangential velocity 
with an amplitude in the range between 5 and 500 km s _1 at 
the distance of 2000AU. As described analytically below, we 
expect no HVSs to be produced at larger tangential veloci- 



Note that the original ejection speed of SDSS 
J090745. 0+024507 should have been higher than its observed 
speed because of its deceleration in the Galactic potential. 
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Figure 1. Ejection speed v e j of the unbound HVS as a function 
of the initial tangential velocity v± of the binary at a distance 
of 2000 AU from SgrA*. The dashed line indicates the analytic 
model of Eqs. J1J and Each point represents the average ejec- 
tion velocity for all HVSs at the given tangential velocity (aver- 
aged over the initial orbital phase within the binary). The aver- 
age value for v c j over all points agrees closely with our analytical 
model. For a = 0.05 AU we get agreement within 16%, and for a 
= 0.10 AU we get agreement within 3%. 



ties. We ran 300 cases for each of the two binary separations, 
with a total of 600 simulations. 



3 ORIGIN OF THE HYPERVELOCITY STAR 

Given a binary system with stars of equal mass m separated 
by a distance a and a SMBH of mass Af > m at a distance b 
from the binary, tidal disruption would occur if b < b t where 
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The distance of closest approach in the initial plunge of the 
binary towards the SMBH can be obtained by angular mo- 
mentum conservation from its initial transverse speed v± at 
its initial distance from the SMBH, d, 

Ld =(^Y ,2 b . (2) 
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The binary will be tidally disrupted if its initial transverse 
speed is lower than some critical value, 
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where a_; 

(m/3Mg), and we have adopted M = 4 x 10 s M Q . For 
v± < fx.crit, one of the stars receives sufficient kinetic en- 
ergy to become unbound, while the second star is kicked 
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into a tighter orbit around the SMBH. The ejection speed 
v c j of the unbound star can be obtained by considering the 
change in its kinetic energy ~ vSv as it acquires a velocity 
shift of order the binary orbital speed Sv ~ yj Gm/a dur- 
ing the disruption process of the binary at a distance ~ b t 
from the SMB H when the binary center-of-mas s speed is 
v ~ y/GM/bt <Hilldll988t lYu fc Tremahiell2003l) . At later 
times, the binary stars separate and move independently 
relative to the SMBH, each with its own orbital energy. For 
v < fx.crit, we therefore expect 

= 1.7 x 10 3 mJ / s 3 alJ /2 km s -1 . (4) 

Figure compares the above approximate model 
(dashed line) with the results from our N-body simulations 
(points). The expected values of w e j and Vx, C rit (dashed line) 
and their dependence on the binary separation a in Eq. J3J, 
are consistent with our numerical results. However, statisti- 
cal variations exist. For v± < ux, cr it the numerical runs show 
variations by up to a factor of ~ 2 around the expected flat 
value of fej in Eq. (QJ . There are also exceptions of escaping 
stars with v± > 200 kms -1 . 

Hills (1988) concluded that the ejection speed of a HVS 
could reach a value of ~ 4000 km s -1 . Although the vast ma- 
jority of our simulated HVSs did not go beyond Hills' limit, 
there were a few exceptions (see Figure Of our 600 runs, 
there were 307 stars that escaped the SMBH, and of those 12 
had velocity 4000 kms -1 . Furthermore, 3 had v e j > 6000 
kms - , of which the largest had velocity v e j = 7073 kms -1 . 
The likelihood of observing such a HVS is remote. As noted, 
HVSs with v e j ^ 4000 kms - rarely occur, and moreover a 
star with velocity 6000 kms -1 would traverse the ~ 200 kpc 
scale of the Milky Way halo \ Wilkinson fc Evanslll99gl) in 
just 30 million years, reducing the likelihood for the observer 
to find it. 



4 FATE OF THE COMPANION STAR 

The orbits of a number of known stars around Sgr A* 
have been studied in detail (see e.g. |Eckart fc Genze1ll997t 
ISchodel et all 120031: iGhez et al.ll2005h . An intriguing ques- 
tion is how did some of these stars obtain their highly ec- 
centric orbits near the central SMBH. 

Figure [3] shows the orbits of the companion stars for 
five HVSs produced by our simulations. All of the orbits 
have very high eccentricities, ranging from e = 0.966 to 
e = 0.999. Similar results are obtained for other orbits, 
not shown in the figure. Previously derived constraints on 
HVS J090745. 0+024507 jFuentes et all2005l) indicated that 
its former companion must have remained bound to SgrA* 
with an eccentricity within the range e = 0.97 to e ~ 1, in 
agr eement with our res ults. 

IGhez et al.l i2005h list stellar orbital parameters for 
seven stars near SgrA*. All but one of the stars have high 
eccentricity, and in particular, SO-16 has a well defined ec- 
centricity of 0.974 ± 0.016. It has a semimajor axis a — 1680 
± 510 AU, and period P = 36 ± 17 years. Our runs produce 
a smaller semimajor axis and period, with a ~ 400 AU and 



P ~ 4 years, but these numbers can be changed by vary- 
ing the stellar masses (see the b ottom right panel of Fig. [3J 
or the initial binary eccentricity. ISchodel et a l. 2003 provide 
the eccentricity of six stars; of particular interest are S14 
(=SO-16) and S8 (=SO-4) with high eccentricities of 0.97 
and 0.98 respectively. The estimated periods are 69 years 
and 342 years for semi major axes o f 3115 AU and 5982 AU, 
respectively. For S0-2. lGIiez et al] i2005l) provided a period 
of 14.53 years and a semimajor axis of 919 AU, with an esti- 
mated eccentricity of 0.8670, which is not wi thin the range of 
our runs but close. iGould fc Quillenl i2003fl suggested that 
the tidal disruption of a massive-star binary could account 
for the orbit of star S0-2. If so, a star of order 100 Mg must 
be the companion. Although possible, such an association is 
unlikely give n the rarity of expected companion stars with 
this mass rsee lKroupal2005h . Stars on an eccentric orbit may 
also be produced through an exchange reaction of a massive 
star wi th a stellar-mass black h ole on a tight orbit around 
SgrA* l|Alexander fc Lividl2004Tl . 

Owing to the flux limit inherent in infrared observa- 
tions of the Galactic center, the observed close-in stars near 
SgrA* are more massive than we assumed in our analy- 
sis. For example, IGhez et al] (I2003T) estimated a mass of ~ 
IOMq for SO-2. An opposite selection effect applies to SDSS 
J090745. 0+024507, because stars with a mass > 10M Q 
would not be observable in the Galactic halo as their life- 
time would be shorter than the duration of their journey. 
The scaling in Eq. Q implies only a modest change in v c j 
for IOMq stars. As an example, the bottom right panel of 
Figure |3] shows the outcome of the disruption process of 
a tight (a = 0.2AU) binary containing stars of IOMq and 
3Mq. The 3Mq star is ejected while the IOMq companion 
remains in a highly eccentric orbit with e = 0.999. The or- 
bital semimajor axis is 700 AU, and the period is 9.3 years. 
A different run with a 0.4AU binary starting at 4000 AU 
with a transverse speed of 5 km s -1 left the IOMq com- 
panion with an orbital semimajor axis of 1430 AU and an 
eccentricity of e = 0.99, closer to the observed parameters 
of S14. 

What is the likelihood for a collision between the stars as 
a result of the kick they acquire from their interaction with 
the SMBH? In our runs, the binary was taken to lie in the 
same plane as its orbit around the SMBH. Assuming that 
the impulsive kick is given by the SMBH towards a random 
direction within this plane, the probability for a collision in 
a case that otherwise would have produced a HVS is four 
time the radius of a star (which is ~ 0.01 AU for a 3Mq star) 
divided by the circumference of a circle with a radius equal 
to the binary separation. For a = 0.05AU and 0.1 AU this 
would imply a collision probability of 12.73% and 6.37%, re- 
spectively. Our runs gave consistent statistical results, with 
a collision fraction among HVS orbits of ~ 7.7 ± 2.1% for 
a = 0.05AU and 5.6 ± 1.9% for a = 0.1AU. The likelihood 
for a collision is expected to be smaller in the more general 
case when the binary lies in a different plane than its orbit 
around the SMBH. Averaging over all random orientations, 
the collision probability is [7r(0.02A£/) 2 / 47r (0-05A[/) 2 ] = 4% 
for a = 0.05AU and 1% for a = 0.1 AU. The small likelihood 
introduces only a minor correction to our earlier statisti- 
cal results (which were obtained by approximating the stars 
as point-like particles). An example of an orbit that leads 
to a collision is shown in the top-right panel of Figure [3] 
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Figure 2. Ejection speed v e j as a function of initial orbital phase within the binary (in degrees). The four cases in (a) refer to a binary 
separation a = 0.05 AU. The cases in (b) refer to a = 0.1 AU. The ejection velocity tends to lie within the range of 2000-4000 km s _1 , 
but occasionally goes beyond it. The ejection speed is sensitive to the orientation of the binary when it gets disrupted, which in turn 
depends on the initial phase as well as the other (internal and bulk) orbital parameters of the binary. The different point types refer to 
different choices for the transverse speed of the binary relative to the SMBH at their initial separation of 2000AU. 
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Figure 3. Orbits of the former companions to HVSs. The SMBH is located at the origin, and the binary system is given an initial 
tangential velocity v± along the negative x-axis at an initial distance of 2000 AU along the positive y-axis. The escape trajectory of the 
HVS is delineated by the dashed line. The binary is often disrupted upon its first revolution around the black hole as seen in the bottom 
right and middle two plots. However, sometimes it takes an extra revolution as seen in the other cases. This is most likely caused by 
the changing phase of the binary at closest approach. Note that the x and y axis have different label scales. The top right panel shows 
the orbit of the binary system until collision. The bottom right panel shows a IOMq star and a 3Mq companion with initial orbital 
separation of a = 0.2 AU. In all other cases both stars were 3Mq and had orbital separation a = 0.05 AU. 
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The two stars would not merge as a result of the collision 
if their relative speed significantly exceeds the escape speed 
from their surface (~ 500 km s _1 ). In the example shown, 
the relative speed of the stars at impact was 666 km s _1 . 
Grazing-incidence collisions, which are more probable than 
head- on collisions, may lead to H VSs which are rapidly spin- 
ning dAlexander fc KumarlkoOll) . 

Finally, we note that for the typical impact parameter 
that leads to the break-up of the binary with a ~ 0.05AU 
by the black hole (Eq. 1), the tidal force on a star is less 
than a few percent of the gravitational force that binds the 
star; however, some rare encounters (with a replaced by the 
stellar radius in Eq. 1) may lead to the tidal disruption of 
the stars. 



5 CONCLUSIONS 

Our N-body simulations indicate that tight binaries with 
separations 0.05-0. 1AU which approach within a distance 
< 10AU from SgrA* could produce HVSs with velocities al- 
most an order of magnitude greater than the observed veloc- 
ity of HVS SDSS J090745.0+024507. The orientation of the 
binary at closest approach plays an important role in deter- 
mining whether the binary is tidally disrupted by the SMBH 
and what is the eventual ejection velocity of the unbound 
star in that case (see Fig. |5J. The phase sensitivity origi- 
nates from the fact that the crossing time of the distance at 
which tidal disruption occurs, (GM/6 3 )" 1 ' 2 is shorter than 
the binary orbital time \/2*pi* (Gm/a 3 ) -1 ^ 2 . Numerically, 
we have not been able to identify a simple trend for the 
phase angle at the breakup radius that would appear more 
organized than the results in Figure [5] 

The former companion star to a HVS is typically kicked 
into a highly eccentric orbit with an eccentricity ^ 0.966 (see 
Fig- El - The resulting eccentricity is similar to that observed 
for a number of stars near the Galactic center such as S8 
(e = 0.98) and S14 (e = 0.97), suggesting a possible binary 
origin for these stars. 

Surveys of HVSs in the halo of the Milky Way galaxy se- 
lect moderate-mass (m < 5Mq) stars with lifetimes longer 
than their travel times (> 10 8 years) while infrared sur- 
veys of the vicinity of SgrA* select for bright massive stars 
(m > WMq) with short lifetimes (< 10 7 years). The findings 
of existing surveys suggest that both populations of stars co- 
exist. Future extensions of this work may examine a larger 
set of binaries with various stellar masses and internal or- 
bital eccentricity. 
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